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The dielectron widths of T(nS) (n — I, ... ,7) and vector decay constants are calculated using the 
' Relativistic String Hamiltonian with a universal interaction. For T(nS) (n = 1,2,3) the dielectron 

widths and their ratios are obtained in full agreement with the latest CLEO data. For T( 10580) 
f"^ . and T(11020) a good agreement with experiment is reached only if the 45* — 3D mixing (with a 

£\J ' mixing angle 6 = 27° ±4°) and 65* — 5D mixing (with 6 = 40° ±5°) are taken into account. The 

^ , possibility to observe higher "mixed D-wave" resonances, Y(n 3 Di) with n = 3,4, 5 is discussed. In 

■ particular, T(~ 11120), originating from the pure 5 3 Di state, can acquire a rather large dielectron 

X^/^ ' width, ~ 130 eV, so that this resonance may become manifest in the e + e~ experiments. On the 

m 



contrary, the widths of pure D-wave states are very small, r ee (n 3 Di) < 2 eV. 



I ■ 1. INTRODUCTION 

The spectrum of bottomonium is very rich with a large number of the levels below the BB threshold. Among them 
three well established T(n 3 Si) (n — 1,2,3) mesons [l| and the l 3 D 2 state, discovered in 01. There exist numerous 
studies of low-lying bottomonium levels, where different QCD motivated models are used 0,13) IE IE S B @ • However, 
a number of these levels have not been observed yet. 
• Observation of the D-wave states which lie below the open beauty threshold, is a difficult experimental task, as 
t^J- | demonstrated in the CLEO experiment , where to discover the 1 3 £>2 level, four-photon cascade measurements in 
\& ■ the T(3S) radiative decays have been performed. In particular, neither the l 3 Di state nor the members of the 2D 
multiplct, for which potential model calculations give masses around 10.45 GeV [1, HE], i.e., below the BB threshold, 
(f-) . are as yet observed. One of the reasons for that is the very small dielectron widths of pure n 3 Di states (for any 
n): here and in (Io| their values ~ 1 eV are obtained. For that reason an observation of a pure D-wave vector state 
O s \ directly in e + e~ experiments seems to be impossible at the present stage. 

However, the situation may change for the D-wave vector states which lie above the open beauty threshold. For 
these bottomonium states the dielectron widths may become larger, as happens in the charmonium family, where the 
k~j ' experimental dielectron width of •0(3770) (which is only 30 MeV above the DD threshold) is already ten times larger 
than for a pure 1 3 Di state [TJ E3- Moreover, the width of the 2 3 Di resonance "0(4160) is almost equal to that of 
-0(4040), which therefore cannot be considered as a pure 3 3 Si state. Such an increase of the dielectron width of a 
D-wave vector state and at the same time a decrease of the width of an S'-wave state can occur if a rather large S — D 
mixing between both states takes place [III EE EH • 

A theoretical study of the <S* — D mixing between vector states is more simple in bottomonium than in the charmo- 
nium family, since the experimental dielectron widths are now measured for the six states X(nS) (n = 1, 6) [J. It 
is also essential that in the recent CLEO experiments the dielectron widths of low-lying levels, T(nS) (n — 1,2,3), 
and their ratios were measured with great accuracy [lij . These three levels can indeed be considered as pure S'-wave 
states, because for them the S — D mixing is possible only via tensor forces, which give very small mixing angle 
(see the Appendix). Then these pure S'-wave states can be studied in the single-channel approach (SCA). Here in 
particular, we use the well-developed relativistic string Hamiltonian (RSH) [15|]. Moreover, just for these levels a 
comparison of experimental and calculated dielectron widths and their ratios can be considered as an important test 
of the theoretical approach and also of the calculated wave functions (w.f.) at the origin. 
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There are not many theoretical studies of higher bottomonium states [T(J HH, Gjl EH ■ Strictly speaking, for this 
task one needs to solve a many-channel problem, knowing the interactions within all channels and between them. 
Unfortunately, this program is not realized now, although some important steps in this direction have been done 
recently in [191 ] , where a theory of the interactions between the channels due to a strong coupling to virtual (open) 
BB (B S B S ) channels was developed. 

In this paper, using the RSH for the low-lying states Y(nS)(n = 1,2,3), we obtain a good description of the 
dielectron widths and their ratios. After that we apply the same approach to higher bottomonium states (above the 
BB threshold), where the accuracy of our calculations is becoming worse: in particular, within the SCA one cannot 
calculate the mass shift of a higher resonance, which can occur owing to coupling to open channel(s). However, for 
the dielectron widths it is most important to define the w.f. at the origin with a good accuracy. 

There exist several arguments in favor of the validity of the the SCA for higher bottomonium states: 

1. First, in charmonium this approximation gives the masses and the dielectron widths of ^(3770), ^(4040), 
?/>(4160), and -0(4415) with a good accuracy, providing a self-consistent description of their dielectron widths 

dam. 

2. Secondly, an open channel (e.g. BB) can be considered as a particular case of a four-quark system QQqq, and, 
as shown in [2(|, this open channel cannot significantly change the w.f. at the origin of heavy quarkonia, because 
the magnitude of a four-quark w.f. at the origin is two orders smaller than for a QQ meson. 

3. Finally, the masses of higher bottomonium states, calculated in SCA, appear to be rather close to the experi- 
mental values, giving a difference between them equal to at most 50 ± 15 MeV. 

For pure n 3 D\ bottomonium states (any n) our calculations give very small dielectron widths, < 2 eV, while for 
higher D-wave vector states their dielectron widths can increase owing to S — D mixing through an open channel(s). 
Here an important point is that the mass of a higher n 3 Di state (n > 3) appears to be only 40-50 MeV larger than 
that of the T((n+ 1)5) state [|, thus increasing the probability of the S — D mixing between these resonances. 
We will define the mixing angle here in a phenomenological way, as it was done in charmonium [ll| , fl3j . 

Owing to the S — D mixing the dielectron widths of mixed D-wave bottomonium resonances appear to increase by 
two orders of magnitude and reach ~ 100 ± 30 eV, while the dielectron widths of the initially pure n 3 Si resonances 
(n = 4, 5, 6) decrease. We also calculate the vector decay constants in bottomonium and briefly discuss the possibility 
to observe "mixed D-wave" resonances in e + e~ experiments. 



2. MASS SPECTRUM 



The spectrum and the w.f. of the bottomonium vector states (L = 0, 2) are calculated with the help of the RSH 
and a universal static potential from [2l[ . This Hamiltonian has been successfully applied to light [HJ and heavy-light 
mesons [23|, [24j, and also to heavy quarkonia @, HH, Ull- In bottomonium the spin-averaged masses M(nL) of the nL 
multiplets are determined by a simpler mass formula than for other mesons, because it does not contain the self-energy 
and the string contributions which in bottomonium are very small, < 1 MeV, and can be neglected 0- As a result, 
the mass M(nL) just coincides with the eigenvalue (e.v.) of the spinless Salpeter equation (SSE) or with the e.v. of 
the einbein equation, derived in the so-called einbein approximation (EA) [231 ] - 

Here we use the EA, because the nS-wetve functions, defined by the EA equation, has an important advantage as 
compared to the solutions of SSE: they are finite near the origin, while the nS-w.f. of the SSE diverge for any n and 
have to be regularized (e.g. as in [4|), introducing unknown parameters. At the same time the difference between the 
EA and SSE masses is small, < 15 MeV (the e.v. in the SSE are always smaller than the e.v. in the EA) and can 
be included in the theoretical error. In Table [T] the centroid masses M(nL) are given both for the SSE and the EA, 
calculated for the same set of input parameters. 

We do not discuss here the singlet ground state 77b, recently discovered by BaBar (2?j], because the hyperfine (HF) 
interaction introduces extra, not well-known parameters, while our goal here is to describe the bottomonium data, 
not introducing extra parameters and using a universal potential which contains only fundamental parameters - the 
QCD constant A and the string tension. 

The masses of T(n 3 Si) are very close to the centroid masses M(nS) (with the exception of the ground state), 
because the HF splittings are small [26|: for higher radial excitations the difference between M{n 3 S\) and the 
centroid mass is < 4 MeV for n > 3. Moreover, for the .D-wave multiplets the fine-structure splittings are small Q 
and therefore the calculated centroid masses M(nD) coincide with M{n 3 Di) within the theoretical error. 

The RSH is defined by the expression from [15|, [22j : 

Ho = ^ + < +uj + VB{r) _ (1) 
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Here nib is the pole mass of the b quark, for which the value nib = 4.832 GeV is used. This number corresponds to 
the current mass mj, = 4.235 GeV, which coincides with the conventional current &-quark mass, equal to 4.20 ± 0.07 
GeV The static potential V B (r), defined below in Eq.©, contains the symbol B, which shows that this potential 
was derived in background perturbation theory [l5| . 

In {!]) the variable lu can be defined in two ways: If the extremum condition is put on the Hamiltonian Ho, lo 
is equal to the kinetic energy operator, lu = -\/p 2 + m\ . Substituting this operator lu into Hq, one arrives at the 
well-known SSE. However, the S-w&ve w.f. of the SSE diverge near the origin and for their definition one needs to 
use a regularization procedure, in this way introducing several additional parameters. 

Instead we prefer to use the EA, where the variable lu is determined from another extremum condition, put on 
the e.v. M(nL). Then Lu(nL) is not an operator anymore, but is equal to the matrix element (m.e.) of the kinetic 
energy operator and plays the role of a dynamical (constituent) quark mass. This constituent mass Lu(nL) grows with 
increasing quantum numbers and this fact appears to be very important for light and heavy- light mesons [23, [2J] , 
while in bottomonium the difference between the dynamical mass w n L and the pole mass nib is not large, changing 
from ~ 170 MeV for the IS ground state up to ~ 300 MeV for higher states, like 6S. 

In the framework of the EA the w.f. of heavy-light mesons have been calculated and successfully applied to determine 
the pseudoscalar decay constants of the D, D s , B, and B s mesons, giving a good agreement with experiment [241 ]. 

It is of interest to notice that in bottomonium the masses, calculated in EA and SSE, and also in the nonrelativistic 
(NR) case (where Lu^n(nL) = mb for all states), do not differ much, even for higher states: such mass differences are 
< 40 MeV (see below and Tables [T] and ■ Still, the w.f. at the origin, calculated in EA, takes into account the 
relativistic corrections and gives rise to a better agreement with the experimental dielectron widths than in the NR 
approach. 

In EA the masses M(nL) are defined by the mass formula: 

o 

m 

M(nL) = uj nL + ^-+ E nL {tu nL ), (2) 
where Lu(nL) and the e.v. E tl l have to be defined solving two self-consistent equations (T3, l23l], namely 



PL 



V B {r) 



Vnh{r) = E nL Lp nL {r), (3) 



and the equation 



duj n L 

In ([3]) we use for all mesons the universal static potential Vb (r) from 0, [H| : 



2 2 dE nL 

u nL =m b -^ ■ (4) 



V B {r) =a{r)r - , 5 

3 r 



with the following set of the parameters: 

m b = 4.832 GeV, A B (n f = 5) = 0.335 GeV, 



(6) 



M B = 0.95 GeV, a = 0.178 GeV 

The QCD (vector) constant A B , which determines the vector coupling constant a B {r) (see Eq. ([7]) below), depends 
on the number of flavors and can be expressed via the QCD constant in the MS regularization scheme; the connection 
between both constants has been established in [U,|2^|. In particular, the two-loop constant A#(n/ = 5) = 335 MeV 
in Eq. ([6]) corresponds to the two- loop A^g = 244 MeV, since they are related as As(n/ = 5) = 1.3656 A^jg [28l |. 
However, one cannot exclude that for low-lying bottomonium levels, like IS, IP, and ID, the choice nt = 4, equal to 
the number of active flavors, might be preferable, giving for their masses a better agreement with experiment. Here 
for simplicity we take nj = 5 for all states, because we are mostly interested in higher states, above the open beauty 
threshold. The constant Co occurs in the expression for the variable string tension a(r) given by Eq. f8]). 

The vector coupling in the coordinate space is defined via the strong coupling in the momentum space ct B (q) 



2 f sm(<2T 
ot B \r) = - / dq a B {q), 

7T J q 



II 



a B(q) = i-r 1_ ^2-i — ( 7 ) 



(3ot B V Po t B 
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TABLE 1: The spin-averaged masses M(nL) (MeV) of low- lying multiplets, calculated in nonrelativistic (NR) case, for spinless 
Salpeter equation (SSE), and in einbein approximation (EA). In all cases the parameters of Vb(t) are taken from Eq. ((6]) 



orate 


MR 
IN ±1 




F A 


Pvn Ml 


ID 


y40y 


o/i p;q 
y4oo 


y4oz 


O/I fin QO-l-O OR /l3c \ 


IP 


9894 


9884 


9888 


9900.1±0.6 


25 


10028 


10010 


10021 


10023.3±0.3 (2 3 Si) 


ID 


10153 


10144 


10146 


10161. 1±1. 7 (l 2 Di) 


2P 


10270 


10256 


10261 


10260.0±0.6 


IF 


10355 


10345 


10347 




35 


10379 


10356 


10369 


10355.2±0.5 (3 3 5i) 


2D 


10460 


10446 


10450 




3P 


10562 


10541 


10551 





witht B = ^£l). 

B 

The solutions of Eq. ([3]) are calculated here considering two types of confining potential in ([5]) : one with the string 
tension equal to a constant, cto = 0.178 GeV 2 , and the other with the string tension <j(r) dependent on the QQ 
separation r. Such a dependence of the string tension on r appears if the creation of virtual light qq pairs is taken into 
account, causing a flattening of the confining potential at large distances, > 1.0 fm. This effect may become important 
for bottomonium states with R(nL) > 1.0 fm, giving a decrease of the masses (e.v.). The explicit expression of a(r) 
is taken here from [29| . where it was deduced from the analysis of radial Regge trajectories of light mesons: 

a(r) =a (l- 7 /M) (8) 

with the parameters taken from [2^]: 7 = 0.40, f(r — > 0) = 0, f(r — > 00) = 1.0. 

TABLE 2: The spin-averaged masses M(nL) (MeV) of higher bottomonium states in the NR case, for the SSE, and in einbein 
approximation (EA) for the potential Vs(r) ((5| 



State 


NR 


SSE 


EA 


Exp. [1] 


2F 


10623 


10607 


10613 




45 


10657 


10630 


10645 


10579.4±1.2 (4 3 5i) 


■W 


10717 


10698 


10705 




iP 


10808 


10783 


10795 




3F 


10857 


10835 


10844 




55 


10894 


10862 


10880 


10865±8 (T(10860)) 


4£> 


10942 


10916 


10928 




5P 


11024 


10998 


11009 




65 


11100 


11067 


11084 


11019±8 (T(11020)) 


5D 


11139 


11109 


11123 




IS 


11278 


11240 


11262 




(W 


11310 


11270 


11295 





In Tables [T] and [5] the masses M(nL) are given only for the flattening potential |S]): for low-lying levels they coincide 
with the masses calculated using the linear potential (with a = const = ctq) within < 2 MeV. For higher states and 
using the flattening potential, the masses (e.v.) are smaller by ~ 10 — 60 MeV (see the numbers in Table [3]). In 
particular, the mass difference is only 12 MeV for the 45 and 3D states and already 40 MeV for the 65* and 51? states, 
reaching 64 MeV for the 75* state. It is evident that for a flattening potential the masses M(nS) (any n) are closer 
to the experimental values. 

For a comparison in Tables Q] and [2] the masses M^B,{nL), calculated in the NR approximation (where 
w(nL)=const=rnb for all states) are also given for the same static potential. These masses M^(nL) are always 
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TABLE 3: The masses of higher bottomonium states (MeV) for the static potential Eq. <(Sj with the parameters ((6JI and two 
confining potentials: linear with ao = 0.178 GeV 2 and the flattening potential Eq. Q 



State 


2D 


4S 


3D 


4P 


5S 


AD 


6S 


5D 


7S 


Linear 


10456 


10656 


10717 


10812 


10901 


10950 


11122 


11163 


11326 


Flatt. 


10450 


10645 


10705 


10795 


10880 


10928 


11084 


11123 


11262 



10—20 MeV larger than in the EA, but in its turn the EA masses are 10—20 MeV larger that the e.v. of the SSE. 
Such a small difference between the EA and SSE masses is taken into account here by including it in the theoretical 
error. 

For our further analysis it is also important that due to the flattening effect the w.f. at the origin (for higher states) 
are becoming significantly smaller than for the linear potential, providing a better agreement with the experimental 
dielectron widths. 

Sometimes the point of view is taken that in bottomonium the nonpcrturbative effects (determined by the confining 
potential) play an insignificant role for low-lying levels. To clarify this point we have compared two m.e. for a given 
nS state: of the confining (nonperturbative) potential < a(r)r > and of the gluon-exchange (GE) (or "perturbative" ) 
potential, respectively, introducing their ratio rj(nS): 

( n\ < Cr ( r ) r >nS fn s 

v(nS) = —— — . (9) 

< l%E( r )l >»s 

The results of our calculations are presented in Table [U 



TABLE 4: The ratios v(nS) 



State 


IS 


2S 


3S 


45 


5S 


6S 


v(nS) 


0.24 


0.93 


1.80 


2.78 


3.87 


5.12 



The values of rj(nS) from Table [4] show that only for the IS" ground state the nonperturbative contribution is rather 
small, equal to 24%, while already for the 25 state both contributions are equally important. For higher nS states 
the nonperturbative contribution dominates, being ~ (n— 1) times larger than the perturbative one. For that reason 
the GE potential can even be considered as a perturbation for higher resonances. 

We estimate the accuracy of our calculations to be equal to 15 MeV. The calculated masses weakly depend on the 
admissible variations of the parameters taken (the same accuracy was obtained in studies of heavy-light mesons [24[ 
and the charmonium family [12j|). Still, for higher resonances the accuracy of the calculated masses is worse, since the 
influence of open channel(s) is not taken into account. Here we can only estimate possible hadronic (decay channel) 
shifts, while comparing calculated and experimental masses: for T(10580) and T(11020) a downward shift ~ 50 ± 15 
MeV is expected, while the mass M(5S), calculated in single-channel approximation, is close to the experimental 
mass of T(10860) (see Tabled]). 

Up to now, many bottomonium states, even those which lie below the BB threshold, have not yet been discovered, 
among them the ID multiplet (two states), the 2D and IF multiplets, and maybe, the 3P multiplet, for which the 
centroid mass M(3P) = 10550(15)MeV 1 , very close to the threshold, is predicted (see Table [T]). The observation of 
these "missing" levels would be very important for the theory. 

For further analysis it is also important that the differences between the masses of the (n + 1)5 and nD states 
(n > 3) are small, decreasing for larger n: their values are equal to 60, 48, and 39 MeV for n = 3, 4, and 5, respectively 
(see Table©. 

The w.f. of the nS and nD states are given in the Appendix, together with m.e. like u>(nL), (p 2 ), and those 
which are needed to determine the dielectron widths and vector decay constants. Also we estimate the relativistic 
corrections, calculating the velocities v 2 /c 2 for different states: their values do not change much, from 0.07 for T(15) 
up to 0.11 for T(65) (see the Appendix). These numbers illustrate the accuracy of the NR approximation. 



Here and below in the brackets we give a theoretical uncertainty 
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In conclusion we would like to stress two points again: first, in bottomonium the centroid masses M(nL) coincide 
with the e.v. of the dynamical equation; secondly, the nonperturbative dynamics dominates for all T(nS) with n > 2. 



3. DIELECTRON WIDTHS 

The dielectron widths are defined here with the help of the van Royen— Weisskopf formula [3(| and taking into 
account the QCD radiative corrections [3l|. The widths T ee (nS) and T ee {nD) can also be expressed through the 
vector decay constants fy, for which explicit expressions were derived in the framework of the field correlator method 
in (24J. For the 5- wave states we have 

r ee (n 3 S!) = -^^fy(nS)py = -^\R nS (0)\^ nS (3v , (10) 
and a similar expression is valid for the -D-wave vector states: 

Airp^pp" Ap^Pv^ 

r ee (n 3 A) = ^^fy(nD)0 v = W-\R nD (Q)\^ nD (3 v , (11) 

oM nD M nD 

if the Z?-wave w.f. at the origin is defined according to the expression Eq. (fl"4)) below, which was derived in [32| . 
In Eqs. (fTU|) and (|11[) the QCD one-loop perturbative corrections enter via the factor [3y (31}: 

16 

0v = l-—a.(M v ). (12) 

However, one cannot exclude that higher order perturbative corrections may not be small and therefore, strictly 
speaking the factor /3y, as well as a s (My) in Eq. (TT2"|) . has to be considered as an effective constant. Nevertheless 
this factor cannot be used as an arbitrary parameter. In different approaches its value typically varies in the range 

0. 75 ± 0.05 [1, [H, [Hj], which corresponds to an effective coupling a s '{My) = 0.14 ± 0.04. (About the choice of the 
renormalization scale, taken here equal to the mass of a vector bb meson My, see the discussion in [13]). Here we 
will neglect in the scale the difference between the mass values for higher states, since all of them lie in narrow range, 
10.6-11.1 GeV. 

As a first step we analyse here the dielectron widths of low-lying levels r ee (T(n5)) (n — 1,2,3) and their ratios 
r(m/n), because these do not depend on the factor fty. As a second step, the values of (3y are extracted from 
the magnitudes of the dielectron widths, which are now known with great accuracy owing to the CLEO data [ijj]. 
Surprisingly, just the same value (3y — 0.80 ± 0.01 is extracted from our fits to three dielectron widths T ee (nS) (n = 

1, 2, 3). This value of (3y corresponds to o: s (~ 10.6 GeV) = 0.12 ± 0.01, which appears to be ~ 15% smaller than the 
strong coupling a s (10.330 GeV) = 0.142 ± 0.056, recently extracted from the CLEO data on the total cross sections 
in e + e" annihilation [35| . 

It is reasonable to assume that such a difference may occur due to second and third order perturbative corrections, 
which were taken into account in the CLEO analysis, while second and higher-order perturbative corrections to the 
dielectron widths are not calculated yet. Taking the central value from [35|], a s (10.330 GeV) = 0.142, one obtains 
/3y = 0.76, which is only 5% smaller than our number (3y(My) — 0.80. From this comparison one can estimate that 
in bottomonium the contribution from unknown higher order corrections to the dielectron width is positive and small, 
< 10%. 

In theoretical studies of the dielectron widths and vector decay constants QCD radiative corrections are often 
neglected, i.e., 0y = 1.0 is taken (3(| H3, H^|, while in our analysis only with {3y — 0.80(1) a good description of 
the dielectron widths is achieved. On the contrary, in [ItJ a significantly smaller number, (3y = 0.46, is exploited. 
Probably, such a small value of /3y (or large strong coupling) has been used in [TtJ in order to suppress the large 
values of the w.f. at the origin for low-lying states, obtained in their model. 

Thus we start with the ratios of the dielectron widths for the n 3 Si states (n = 1,2,3): 

r(m/n) - T ^ mS ^ (M(nSr)(R mS (0)f 
r [m/ n) - r ^ {nS) - {M{mS) 2 ){RnSm 2> 

which are fully determined by the w.f. at the origin (the masses are known from experiment). Taking the w.f. at the 
origin calculated here, from the Appendix, one arrives at the values of r(m/n) given in Table O 

Both the calculated and the experimental ratios agree with each other with an accuracy of < 3% and this result 
can be considered as a good test of our approach, as well as of the w.f. at the origin calculated here. 
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TABLE 5: The ratios of the dielectron widths r(m/n) for low- lying n 3 Si states 





r(2/l) 


r(3/l) 


r(3/2) 


Theory 


0.465 


0.339 


0.728 


Exp. [14] 


0.457 ±0.008 


0.329 ± 0.006 


0.720 ±0.016 



TABLE 6: The dielectron widths of n 3 Si (n = 1, 2, 3) and n 3 Di (n = 1, 2) keV with (3 V = 0.80 



Widths 


Theory 


Exp. [14] 


r ee (is) 


1.320 


1.354±0.024 


r ee (2S) 


0.614 


0.619±0.014 


r ee (35*) 


0.447 


0.446±0.011 


r„(iD) 
r ee (2D) 


0.614xl0" 3 
1.103X1Q -3 





Next we calculate the absolute values of the dielectron widths, which allow to extract the QCD factor f3y. From three 
dielectron widths T ee (nS) (n — 1, 2, 3) the same value By — 0.80(1) has been extracted. Later everywhere f3y = 0.80 
is used, for which the dielectron widths of low- lying and higher states are given in Tables [B] and [71 respectively. 

For low-lying levels the dielectron widths (with (3y — 0.80) agree with the experimental numbers within 3% accuracy 
(see Table©. 

The dielectron widths calculated here are compared with other theoretical predictions [It], in Table [Jj in [ItJ 
rather small dielectron widths are obtained, mostly due to the small By = 0.46 taken there. This value is 70% smaller, 
i.e., the QCD radiative corrections are larger, than in our case. In [l8|, as well as in our calculations, for T(10580) the 
dielectron width is larger than in experiment, while for the ground state their dielectron width is three times smaller 
than in our calculations and in experiment. 



TABLE 7: The dielectron widths F ee (nS) (keV) of pure S- wave states 



State 


IS 


2S 


3S 


4S 


5S 


6S 


GVGV a [17] 


1.01 


0.35 


0.25 


0.18 


0.14 




CO [18] 


0.426 


0.356 


0.335 


0.311 






This paper 


1.320 


0.614 


0.447 


0.372 


0.314 


0.270 


Exp. [1] 


1.340 ±0.018 


0.612 ±0.011 


0.443 ± 0.008 


0.272 ±0.029 


0.31 ±0.07 


0.13 ±0.03 



a The numbers given are taken from second paper in [17|I 

In conclusion of this section we would like to stress again that: 

1. The value By = 0.80(1) should be considered as an effective constant, which implicitly takes into account the 
contributions from higher perturbative corrections. We expect that higher-order perturbative corrections are 
positive and rather small, < 10%). In the absence of higher corrections, the effective coupling, a s (10.6 GeV) ~ 
0.12(1) taken here, appears to be slightly smaller than the strong coupling, extracted from the analysis of the 
cross sections of e + e~ — > hadrons [35j. 

2. The calculated ratios of the dielectron widths (for low-lying levels), which are independent of the unknown QCD 
factor By, agree with experimental ratios with an accuracy better than 3%. Therefore one can expect that in 
our approach the w.f. at the origin (for low-lying levels) are calculated with a good accuracy. 

3. The dielectron widths calculated here (with By — 0.80) agree with experiment with an accuracy better than 
5%. 
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4. THE S-D MIXING BETWEEN THE (n + 1) 3 Si AND n 3 L>i BOTTOMONIUM STATES 

In contrast to the case of the low-lying levels, the calculated dielectron widths of pure nS vector states with n = 4 
and n = 6 exceed the experimental values: for the AS and 6S states they are 25% and two times larger than the 
experimental widths of T(10580) and T(11020), respectively. Such a suppression of the dielectron widths occurs if 
one or more channels are open. Some reasons for that have been discussed in [l6| . where it was that in particular in 
the Cornell coupled-channel model [3!| the dielectron widths of higher charmonium states are not suppressed. 

Here as in jl 31 ] - we assume that an open channel cannot significantly affect the w.f. at the origin calculated in 
closed-channel approximation. This assumption is based on the study of a four-quark system in (20| . where the 
calculated w.f at the origin of a four-quark system, like QQqq, appears to be about two orders smaller than that of 
a heavy meson QQ. We expect this statement also to be true of a continuum w.f. at the origin (the w.f. of an open 
channel), which can be considered as a particular case of a four-quark system (this does not exclude that a continuum 
channel can strongly affect the QQ w.f. at large distances). Thus it is assumed that a suppression of the dielectron 
widths of higher states occurs due to the S — D mixing between the [n + 1)S and nD vector states, which happen 
to have close values of their masses. We also show in the Appendix that in bottomonium the S — D mixing due to 
tensor forces appears to be very small, giving a mixing angle Bj- < 1°. 

For the D-w&ve states their w. f. at the origin is defined here as in (32J: 

e> fm _ 5i%(0) (u) 
RD{0) - (M) 

and for the mixed states their physical (mixed) w.f. are given by 

i?phy S s(0) = cos0ite(O) - sin0 R D {0), (15) 
flphys d(0) = sin 9R S (0) + cos 9 R D {0). (16) 

The w.f. at the origin of pure S- and D-w&ve states and the derivatives (0) are given in the Appendix together 
with other m.e. which are needed to calculate the physical w.f. at the origin (see Table [T4"|) . As seen from Table [T2l 
the w.f. i?„r)(0) are small and therefore the dielectron widths of pure n 3 Di states appear to be very small, < 2 eV. 
Their values are given in Table [8] Notice that our widths are ~ 10 times smaller than those in [36]). On the contrary, 
the dielectron width of the 4 3 Si resonance is 25% larger. 

To obtain agreement with the experimental value, r ee (10580) = 0.273 ± 0.022 keV, we take into account the 45 — 
3D mixing and determine the mixing angle, 9 = 27° ± 4° from this fit. Thus T(10580) cannot be considered a pure 
4S vector state, it is mixed with the initially pure 3 3 Z?i state. This second "mixed" state will be denoted here as 
T(~ 10700), it acquires the dielectron width r ee (T(10700)) = 0.095 keV, which is 60 times larger than the width of 
the pure 3 3 Di state. 

The dielectron widths are given in Table [8] in two cases: without the S — D mixing [6 — 0) and for mixed states, 
taking the same mixing angle 9 = 27° for all higher states. 

An interpretation of the experimental width of T (10860) cannot be done in an unambiguous way: the calculated 
r ee (5 3 5i) of a pure 55* state (9 = 0) just coincides with the central value of the experimental r ee (T(10860)) = 
0.31 ± 0.07 keV. It could mean that for some unknown reason the 5S and AD vector states are not mixed. However, 
there exists another possibility, because the width r ee (10860) has a rather large experimental error. In particular, for 
the mixing angle 9 — 27° one obtains r ee (T(10860)) = 0.23 keV, which just coincides with the lower bound of the 
experimental value. To decide which of the two possibilities is realized, more precise measurements of r ee (T(10860)) 
are needed. 

An interesting opportunity can be realized for the originally pure 5 3 Di resonance. The experimental dielectron 
width of T(11020) is very small, r ee (11020) = 0.13±0.03 keV, i.e., it is two times smaller than the number calculated 
here without the 6,!? — 5D mixing (i.e., 9 — 0). Even for the mixing angle 9 = 27°, the theoretical value is still 26% 
larger compared to the experimental one (see Table [8]). To fit the experimental number a rather large mixing angle, 
40° ± 5°, has to be taken. For a such a large angle the dielectron widths of both resonances, T(5D) (with mass 
~ 11120 MeV) and T(11020), appear to be almost equal: 

( r ee (T(11020)) = 0.137 ±0.025 keV, , . 

[ r ee (f (5Z>)) = 0.135 ±0.025 keV. ^ ' 

It is of interest to notice that this large angle is close to the value of the mixing angle 9 = 35°, which has been 
extracted in [l^] and to fit the dielectron widths in the charmonium family: ip(A0A0), t/>(4160), and ^(4415). 



9 



TABLE 8: The dielectron widths (keV) for the mixing angles 6 = and 9 = 27° (f3 v = 0.80) 



Widths 


Theory 




Exp. [14] 




no mixing with mixing 




T ee (45') 


0.372 


0.273 


0.272±0.029 


r ee (3D) 


1.435X10" 3 


0.095 




r ee (5S) 


0.314 


0.230 


0.31±0.07 


r ee (4£>) 


1.697xl0~ 3 


0.084 




r ee (65) 


0.270 


0.196 


0.13±0.03 


r ee (5D) 


1.878xl0" 3 


0.075 





TABLE 9: The decay constants fv(nS) (MeV) 



State 


IS 


2S 


35 


45 


5S 


6S 


Pv = 1.0 [36] 


498(20) 


366(27) 


304(27) 


259(22) 


228(16) 




This paper, fiv = 0.80 


794 


557 


483 


383 


355 


331 


jv(exp) for (5v = 0.80 [1] 


798 ±6 


556 ±6 


481 ± 5 


381 ± 19 


413 ± 45 


268 ± 30 


/v(exp) for 0V = 1.0 [1] 


715 ±5 


497 ±5 


430 ±4 


341 ± 17 


369 ± 40 


240 ± 27 



5. DECAY CONSTANTS IN VECTOR CHANNELS 

The decay constant in vector channel fv{nL) is expressed via the dielectron width in a simple way, as in Eqs. (|10p 
and (fTTj) . Therefore, from the experimental widths the "experimental" decay constants can be easily obtained. Still 
an uncertainty is left, coming from the theoretical error of about 10% in the QCD factor (3y. Also in many papers 
perturbative one- loop corrections are neglected, i.e., (3y — 1-0 is taken [H, H3, HI] . This makes a comparison with 
other calculations more difficult. In our calculations we take f3y = 0.80, which is slightly more than (3y ~ 0.70, used 
in [| and [H- 

To determine the experimental /V(exp) we take in this section the experimental data from PDG jl| (not the CLEO 
data [3), which are used in most of the cited theoretical papers. These decay constants are given in Table[9l both for 
/V = 1-0 and Pv = 0.80, the difference between them is ~ 10%. The theoretical predictions for fy give significantly 
different numbers [3g, [37], [38| , which shows that the decay constants are rather sensitive to the dynamical parameters 
of the interaction and the model used. For a comparison we give in Tables 191 and [TPl the decay constants fy(nS) and 
fv(nD), calculated here and in Rcf. [36], where the relativistic Bethe— Salpeter method was used. All values needed 
for our calculations are presented in the Appendix. 

For a comparison we also mention here the values of fy(lS) from [3?], HH where the QCD factor f3y = 1.0 was used: 
f v (lS) = 529 MeV in [37| is significantly smaller than the value fvO-S) = 705(27) MeV in [H, which is very close 
to the "experimental" fy(lS) (see Table [9]). On the contrary, in [33| the perturbative corrections have been taken 
into account with f3y ~ 0.66. There the values of fy are not given but the dielectron widths of the nS vector states 
(ri = 1, 2, 3) are in reasonable agreement with experiment. 

For the -D-wave states the w.f at the origin, the second derivatives, and other m.e. determining the vector decay 
constants fy (nD) , are given in the Appendix. The calculated fy(nD) are presented in Table [T0l together with the 
numbers from [36j. Unfortunately, at present there are no experimental data on the dielectron widths for those states. 

The decay constants of pure nD vector states (n = 1,2,3), calculated in our approach, appear to be ~ 10 times 
smaller than those from [361 ]. where the Bethe— Salpeter equation was used, and the reason behind such a large 
discrepancy remains unclear. However, if the 45 — 3-D mixing is taken into account (with 9 — 27°), then the values of 
fv(3D) are close to each other in both approaches. 

6. SUMMARY AND CONCLUSIONS 

In this paper we have calculated the bottomonium spectrum and shown that the masses of the (n + 1)S and nD 
states (for a given n > 2) are close to each other. We also assume here that between these states S — D mixing 
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TABLE 10: The decay constants fv{nD) (MeV) 



State 




ID 


2D 


3D 


4D 


5D 


Pv = 1-0 [36] 




261(21) 


155(11) 


178(10) 






This paper, /3v = 1.0, ( 


) = 0° 


18 


24 


28 


31 


33 


This paper, (3v = 1.0, 6 


i 2T° a 






226 


215 


206 



a In bottomonium the 2S — ID and 35 — 2D states, occuring below the threshold, do not mix via tensor forces, see a discussion in 
Appendix. 



takes place, which allows to describe the dielectron widths of higher states with a good accuracy. There are several 
arguments in favor of such a mixing. 

1. Suppression of the dielectron widths of T(10580) and T(11020). 

2. Similarity with the S — D mixing effects in the charmonium family. 

3. Strong coupling to the BB (B S B S ) channel. This fact has been supported by recent observations of the resonances 
in the processes like e + e~ — > T(nS)n + TT~ (n = 1, 2, 3) and the theoretical analysis in 19]. 

The important question arises whether it is possible to observe the mixed D-wave vector resonances in e + e~ experi- 
ments. Our calculations give M(3D) = 10700(15) MeV (not including a possible hadronic shift) and r ee (T(3Z3)) ~ 95 
eV for the mixing angle 6 — 27°, which is three times smaller than r ee (T(10580)). For such a width an enhancement 
of this resonance in the e + e~ processes might be suppressed, as compared to the peak from the T(10580) resonance. 

The situation remains unclear with the 5S — AD mixing, because the dielectron width of T(10860) contains a rather 
large experimental error and a definite conclusion about the value of the mixing angle, or no mixing at all, cannot be 
drawn. We have considered both cases here, obtaining the mass 10930 ± 15(th) MeV for the AD state. 

It looks more probable to observe the resonance T(5D) (with the mass ~ 11120 MeV), for which the dielectron 
width may be almost equal to that of the conventional T(11020) resonance. However, since the cross sections of e + e~ 
processes depend also on other unknown parameters, like the total width and branching ratio to hadronic channels, 
the possibility to observe a mixed 5D vector resonance may be smaller than for T(11020), even for equal dielectron 
widths. 

Recently new observations in the mass region 10.6-11.0 GeV have been reported [4lll42Tl. The resonance T(10890), 
considered to be identical to T (10860), has been observed by the Belle Collaboration Two resonances in the 

same region, T(10876) and T(10996), have been measured by the BaBar Collaboration [42j, they are supposed to 
be the conventional T(10860) and T(11020). Still there are some differences between the masses and total widths of 
the resonances from ^3, |42j, and the PDG data [l|, so that further analysis of their parameters is needed. Also one 
cannot exclude that an overlap of T(11020) with the still unobserved T(5D) resonance is possible, which can distort 
the shape and other resonance parameters of the conventional T(11020) resonance. 
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Appendix 

The wave functions at the origin and some matrix elements 



We start with the definition of the vector decay constants fv , 
framework of the vacuum correlator method in 1241: 



for which the following expression was derived in the 



fi(nL) = 12 



h/w?(Q)r 

M v (nL) 



3 \Rns(0)\ 
7T M v (nS) 



Hv(nL). 



(A.l) 
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Besides the w.f. at the origin, the vector decay constant contains a relativistic factor £y(nL), which was also defined 
in M: 



W 2 ) 



2uj 2 



(A.2) 



Numerically, £y(nL) is close to unity for all nS and nD vector states. Moreover, for the (n + 1)5 and nD states, 
which we are mostly interested in here, their values coincide (see Table ITT]) . Notice that all £y(nL) differ from unity 
at most by 4%. 



TABLE 11: Relativistic factors £.v(nL) 



State 


IS" 


2S 


3S 


45 


5S 


6S 


ID 


2D 


3D 


4£> 


5D 




0.976 


0.974 


0.972 


0.968 


0.966 


0.963 


0.975 


0.972 


0.968 


0.966 


0.965 



The meaning of the m.e. uj n s and UnD has been discussed in Section 1 and their values are given in Table 1121 
together with the m.e. (p 2 ) and the w.f. at the origin R n s(0) and R n £)(0). Notice that again the values of ui n L 
coincide for (n + 1)S and nD states. Knowledge of the m.e. (p 2 ) and oj(nL) allows to define the velocities v 2 /c 2 for 
different states, which in all cases appear to be rather small, < 0.11. 



TABLE 12: The matrix elements uj nS (GeV), (p 2 > (GeV 2 ), and the w.f. at the origin R n s{0) (GeV 3 ' 2 ) (no mixing) for the 
potential ([S]) 



State 


IS 


2S 


3S 


45 


5S 


6S 


ID 


2D 


3D 


4D 


5D 


w 


5.02 


5.03 


5.05 


5.08 


5.1 


5.12 


5.02 


5.05 


5.08 


5.1 


5.11 


R(0) 


2.529 


1.829 


1.616 


1.506 


1.424 


1.349 


0.0595 


0.081 


0.0956 


0.107 


0.1145 


(P 2 > 


1.873 


1.907 


2.178 


2.459 


2.689 


2.844 


1.821 


2.136 


2.423 


2.661 


2.810 



The w.f. RnD(0) are defined via the derivative R nD (0) according to the definition 
calculations give the following numbers of the derivatives, presented in Table [T3J 



14]) from Ref. [32]. Our 



TABLE 13: The second derivative R^ D {Q) (GeV 7/2 ) 



state ID 


2D 


3D 


4D 


5D 


< D (0) 0.848 


1.167 


1.396 


1.572 


1.692 



An interesting feature of the D-wave w.f. is that in contrast to the w.f. R n s(0), which decrease for larger n, the 
derivatives R nD (0) and R n D(0) grow for higher radial number n (see Tables ITZl and fT3|) . and this effect increases the 
probability of the S — D mixing. 

In bottomonium the w.f. at the origin has been calculated in many papers (l6l. flsL [43l| . Calculated in our approach 
w.f. for the pure nS and nD states, which arc presented in Table [TJl allow to determine the physical w.f. at the 
origin of mixed states for an arbitrary mixing angle. For higher states such physical w.f. at the origin are given in 
Table [TJ 

For T(10580) and T(10860) we take the mixing angle 9 — 27°, which was extracted to fit the dielectron width 
T ee (10580). For T(11020) we take a larger mixing angle, 9 = 40°, to fit its dielectron width. 

Just the numbers from Table [JJ] have been used in Sections 4 and 5 to calculate the dielectron widths and the 
vector decay constants. 

In this paper we do not discuss the pseudoscalar decay constants fp(nS) in bottomonium, although for them the 
expressions similar to (|A.lj) have been obtained in (24[; these constants need a special consideration. 

In conclusion we also calculate the m.e. over the tensor forces, which are defined by the potential Vr(r): 



V T {r) 



ia s (mb) 1 



3u, b 2 



(A.3) 
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TABLE 14: The physical wave functions at the origin i?ph ys (0) (GeV 3/2 ) 



Mixing angle 


9 = 


27° 


6 = 40° 


State 


T(10580) 


T(10860) 


T(11020) 


#phy S (0) 


1.298 


1.220 


0.960 


State 


f (10700) 


f (10930) 


f (11120) 


^Phys(0) 


0.769 


0.742 


0.955 



Notice, that in denominator instead of m\, usually used in different models, in the field correlator method the tensor 
potential (| A.3[) has to contain the squared dynamical mass u(nL) 2 [ll| |44|. The 5 — D mixing between (n + 1)5 and 
nD vector states, due to tensor interaction, is proportional to the nondiagonal m.c. (f~ 3 )(n+i)s.ri.Dj f° r which the 
following numbers were obtained here (see Table [T5]) . 



TABLE 15: The matrix elements <r~ 3 ) (GeV 3 ) 



n 1 


2 


3 


4 5 


(r- 3 )(n+i)s,nD 0.0574 


0.0683 


0.0721 


0.0739 0.0744 



With the help of these numbers one can easily calculate the m.e. over the tensor potential (|A.3[) . considering it as 
a perturbation. These m.e. are indeed very small, < 1 MeV and therefore give rise to small contributions to the w.f. 
at the origin. These small correction to the w.f. at the origin define the mixing angles due to tensor interaction, 9t, 
for different states. We obtain that an admissible mixing angle 9t is small, 6? < 1°. 

The vector decay constants were calculated with the use of Eq. (|A.1|) and presented in Table [16] For higher states 
their values are given in two cases: without mixing and for 9 = 27°. One can see that for the originally pure 45 
and 65 states the constants fv decrease by 15% and 35%, respectively, compared to the case without mixing. For 
originally pure D-wave states their decay constants fy are growing, namely, fy(ZD) and fy(5D) increase eightfold. 



TABLE 16: The decay constants fv (MeV) for pure and mixed states with 8 = 27° 



State 


15 


25 


35 


45 


55 


65 


ID 


2D 


3D 


4D 


5D 


no mixing 


794 


557 


483 


444 


415 


389 


18 


24 


28 


31 


33 


6 = 27° 








383 


355 


331 






226 


215 


206 



In Table [TBI for low-lying mixed states the values of fy(nS) (n=l,2,3) and fy(nD) (n=l,2) are absent, because 
these states are indeed pure states, since for them the 5 — D mixing can occur only via tensor forces, which give very 
small effect 6 T < 1°. 
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